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Abstract 

We examine the role of the elementary isoscalar pion-nucleon scattering amplitude 
in the description of the processes pp — > pp-K^ and pp — > (i7r+ at threshold. We 
argue that the presently used tree level dimension two approximation used in chiral 
perturbation theory is insufficient as input by direct comparison with the vrN scat- 
tering data. We also show that a successful semi-phenomenological boson-exchange 
model does better in the description of these data. The influence of the violation of 
crossing symmetry in the meson-exchange model has to be studied in more detail. 
We stress that further investigations of the process pp dir^ can pave the way to 
a deeper understanding of the pion dominated part of the transition operator. 
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1. The high precision data for the processes pp ppn'^ and pp dir^ in the threshold 
region |T| have spurred a flurry of theoretical investigations. The first data on neutral 
pion production were a big surprise because the experimental cross sections turned out 
to be a factor of five larger than the theoretical predictions. These included the direct 
(Born) graph (fig. la) and the on-shell rescattering (fig. lb), with the pion-nucleon (ttTV) 
T-matrix replaced by the scattering length. Subsequently, it was argued that heavy- 
meson exchanges (fig.lc) might be able to remove this discrepancy P]. On the other hand 
it was shown [|]|]5| that the off-shell behaviour of the full vrA^ T-matrix also enhances 
the cross section considerably. In , the pion-nucleon T-matrix was calculated within a 
semi-phenomenological meson-exchange model, based on an effective Lagrangian fullfill- 
ing certain requirements from chiral symmetry (as discussed below). Recently this result 
has been questioned based on calculations within the framework of tree-level chiral per- 
turbation theory (CHPT) including dimension two operators 0, ^, ^. In these papers, 
CHPT has been used to constrain the long-range pion-exchange contributions. In fact, 
the chiral perturbation theory approach, which abides to all symmetry requirements in- 
cluding crossin^^, is believed to lead to a deeper understanding of the success of the 
meson-exchange picture, as first stressed by Weinberg [0. It is therefore striking that 
the calculations for pp —* ppix^ performed so far lead to a marked difference in the role 
of the so-called rescattering contribution^ which interferes constructively with the direct 
production in the Jiilich model and destructively in the chiral framework, respectively. 
Note that while there is still debate about the actual numerical treatment (co-ordinate 
versus momentum space) and the ensuing size of the rescattering contribution in the chi- 
ral perturbation theory approaches the sign difference to the meson-exchange model 
can be considered a genuine feature. It is exactly this point which we wish to address in 
this letter. We argue that the treatment underlying the isoscalar pion-nucleon scattering 
amplitude and the related transition operator for the process NN —>■ NNir in the chiral 
framework is not yet sufficiently accurate and thus the resulting rescattering contribution 
should be considered an artifact of this approximation. Clearly, this does not mean that 
chiral perturbation theory is invalid but rather that higher order (one loop) effects need to 
be accounted for. This statement can be made very transparent if one considers in detail 
the isoscalar vrA^ scattering amplitude. We stress, however, that in the meson-exchange 
model crossing symmetry is violated and we can not exclude the possibility that this is 
in fact the reason for the constructive interference. To gain a deeper understanding of 
the underlying direct and rescattering contributions, we propose to study in more detail 
the process pp — >■ dn'^ since heavy meson exchanges play no role in that case and one 
thus has a better chance to extract the pionic (chiral) physics which underlies the param- 
eterization of the NN — > NNtt transition operator. We furthermore mention that for 
our investigation mainly focusing on the strength of the S-wave rescattering, the explicit 
inclusion of the A(1232) is not needed. Clearly, for a precise comparison with the data, 
this approximation can not be maintained. Some of the effects of the A are encoded in 
the numerical values of the low-energy constants q. Further dynamical effects related 
to the A only appear at higher orders togther with many other terms. A systematic 
inclusion of such effects goes beyond the scope of this letter. We stress that our intention 



*^Note, however, that at tree level unitarity is violated since all amplitudes are real. 
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is to critically assess the accuracy of the presently available CHPT calculations for pion 
production in pp collisions. 

The important feature of this reaction is the large momentum mismatch between the 
initial and the final nucleon-nucleon state: the initial relative momentum is given by pf^ = 
ttiMtt + whereas close to threshold the final one is compatible with zero. This huge 

difference leads to momentum transfers in the order of pm, at least if one restricts oneself 
to S-wave production. Let us concentrate first on the so-called rescattering contribution 
depicted in fig. |l]b. To get a feeling for the relevant momentum transfer, consider the 
reaction at threshold and neglect the pp final-state interactions. One finds for the four- 
momentum of the exchanged pion, k"^ = —itlMt^ ~ ~7M^ ~ —0.1 GeV^, which is deeply 
space-like. We remark that in single pion electroproduction off protons, one-loop effects 
for such photon virtualities are fairly large |]TT| substantiating our previous statement that 
one has to go beyond tree level in CHPT. 

2. The study of pion-nucleon scattering in chiral perturbation theory has already a 
long history. In the relativistic approach, Gasser et al. [|l^ constructed the full off- 
shell ttN amplitude to one loop in the chiral expansion. Their work focused mainly 
on the determination of the pion-nucleon a-term and the related scattering lengths. 
A systematic comparison with the low energy vrN scattering data was, however, never 
presented. In the framework of heavy baryon CHPT, the chiral corrections to the S-wave 
scattering lengths were worked out by Bernard et al. |T^. In particular, it was shown 
that the isoscalar scattering length a+ can be expressed as (to one loop) 



m J F2 V 8m/ QA-kF^ 

with (^) the charged pion (nucleon) mass, F^^ = 92.4 MeV the pion decay constant and 
qa = 1.26 the nucleon axial-vector coupling constant. The q are low-energy constants 
(LECs) not fixed by chiral symmetry that appear in the dimension two chiral Lagrangian, 

C^^l = NLTrix+) + ic2-^)iv-uf + csu-u + ..)N , (2) 
V 8m / 

where the ellipsis denotes further terms not needed here, is the standard axial-vector, 
which transforms homogeneously under non-linearly realized chiral symmetry, the 
nucleons' four- velocity (f^ = 1), x+ ~ parameterizes the explicit chiral symmetry 
breaking due to the light quark masses (and thus is directly related to the a-term) and 
N is the velocity-projected nucleon isodoublet, = (1+ y^^/2 (for a detailed review, 
see |T^). The q are finite since loop corrections start to contribute only at the next order 
in the chiral expansion. The precise values of these LECs will be discussed below. It 
is important to note that there are large cancelations between the one loop contribution 
of order and the kinematically 1/m corrections at order M^. This makes the actual 
value of a"*" very sensitive to the values of the LECs q. The scattering amplitude in the 
tree approximation to second order in small momenta is readily derived from Eq. (Q) since 
to leading order in pion fields, = —id^j,(j)/F.,t + 0{(lP'). This is the accuracy to which this 
amplitude has been used in chiral perturbation theory approaches to pp —>■ ppn^ in the 
threshold region. In that case, the produced pion is almost at rest, i.e. has a very small 
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three-momentum (denoted q) , whereas the exchanged pion typically carries a momentum 
A; of a few hundred MeV. The explicit form of the half-off shell scattering amplitude 
T^{q,k) is e.g. given in 0. Clearly, the sensitivity to the precise values of the LECs 
observed for the scattering length carries over to the scattering amplitude. Note that this 
isoscalar scattering amplitude gives rise to the so-called rescattering contribution. 

To proceed, we briefly summarize what is known about the LECs q. First, it is 
important to notice that there is an ambiguity in pinning down their values. These 
depend on the order one is working. From a tree level fit to pion-nucleon scattering 
(sub)threshold parameters including the dimension two operators, one obtains [^] 



ci = -0.64 ±0.14 , C2 = 1.78 ±0.10 , C3 = -3.90 ± 0.09 , (3) 

with all numbers given in GeV~^. Note that the sum c"*" = C2 ± C3 — 2ci which enters the 
isoscalar amplitude, is = —0.84. In this determination, the empirical value for a"*" was 
not used simply because it is badly determined and also, it is not possible to describe all 
(sub)threshold data that are sensitive to the q by one simultaneous fit at order (as 
discussed in detail in [|^). We therefore have also performed calculations with = 0.459 
and 0.005 GeV~^ corresponding to the the conservative band of = ±10 ■ lO'^/M^r, 
respectively. Since there are no compelling arguments which of the q should be readjusted, 
we have changed either one of them so as to get these values for a~^. 

To one loop order, one can determine these LECs from a set of seven observables 
which are given by tree graphs including the Cj and finite loop corrections but have no 



contribution from the 24 dimension three LECs [ICl. One finds 



ci = -0.93 ±0.10 , C2 = 3.34 ±0.20 , C3 = -5.29 ± 0.25 . (4) 

The differences between the numbers given in Eq.(0) and Eq.(^) are sizeable showing 
that the one-loop corrections are not small and thus have to be taken into account. We 
remark that in this case = —0.09 and a"*" = —4.7 ■ 10~^/M^, well within the empirical 
band. Note that c"*" is an order of magnitude smaller than any of the LECs individually 
and is well within the uncertainty of the LECs. This observation is at the origin of the 
statement that to understand the empirical value of a~^, one has to know these LECs 
very precisely. Mojzis fl^ has recently performed a complete one-loop calculation of 
the ttN scattering amplitude, including the dimension three counterterms. A fit to the 
known S-, P-, D- and F-wave threshold parameters as well as to the pion-nucleon a-term 
and the Goldberger-Treiman discrepancy allows to pin down all LECs. Remarkably, the 
dimension two LECs of relevance here come out completely consistent with the values 
found in |TB[ . To assess the sensitivity to the LECs, we will use both sets even so we will 



only work to second order in the chiral expansion. This was also done in At this 

point, one might already suspect that this approximation is not sufficient. 

In contrast, the isovector amplitude is dominated by the Weinberg term, i.e. the nirNN 
vertex which stems from the chiral covariant derivative and has it strength given entirely in 
terms of 1/ 4F^. The chiral corrections have also been calculated [|13| [jl8[ , but since they are 
small, they are not important for the following arguments. By forming appropriate linear 
combinations of the isoscalar and isovector amplitudes, one can construct the phase shifts 
Sii and 5*31. They consist of a Born and a non-Born piece and take the form (expanding 
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the partial waves, not the invariant amphtudes, to second order in small momenta) 
Sn = V/-^ - Afl (l - + - + 2Mje, - 2^^(c, + c,) 



2 



4 V 2m 3m 

S3. = v';:;jr]si(i_^)_i_||__L(^^^^^ 



6m 



2 r /\^4 



(5) 



with UJ the pion cms energy. The phase shifts are shown in fig. Q for the Cj as given in 
Eqs.(I^J^) in comparison to the data. In both cases, the prediction deviate already at low 
energies from the data. Due to the approximations involved, in both cases the partial 
waves increase quadratically at higher energies, i.e. one has to include loop effects if one 
wishes to find a good representation of the vrN phase shifts for energies up to 100 MeV or 
higher. We have also calculated these phases at tree level (dimension two) with c"*" adjusted 
so as to give the empirical range of as discussed above. Enhancing e.g. C2 by a factor of 
1.7 in the tree level fit leads to a good description of the phase shifts up to Tcms — 125 MeV. 
Such an enhancement is, however, incompatible with some subthreshold parameters, for 
detailed study see |TB[. We therefore do not consider this a justified way of describing 



elastic pion-nucleon scattering for the energies under consideration. Furthermore, at 
higher energies the unphysical quadratic growth shows up again. We remark that the 
behavior of the phase shifts shown has direct consequences for the transition operator 
relevant for A^A^ NNir as discussed below. 

3. There is a different way of constructing the A^A^ — > NNn transition operator that is 



by means of a meson-exchange model as e.g. the one constructed by Schiitz et al. |T9|. 
We briefly summarize its salient features without going into any kind of detail (see also 
the discussion in ref.|P). The pseudo-potential contains in addition to the nucleon- and 
A-pole and exchange diagrams p and a t-channel-exchanges. However, the latter two do 
not appear as sharp particles but are constructed from the correlated two-pion exchange 
by use of a dispersion integral. The potential is then iterated in a relativistic Lippmann- 
Schwinger equation within the framework of time-ordered perturbation theory to ensure 
unitarity. To get finite results, the vertices are supplemented by meson-nucleon form 
factors. This approach gives a unique prescription to go off shell. The parameters are 
then fitted to the vrA^ phase shifts up to lab energies of about 500 MeV under the constraint 
that the scattering lengths come out to be in the empirical ranges. In addition it turns 
out that the model produces a reasonable value for the vrA^ a-term. As it was stated 
above the value of the isoscalar scattering length is compatible with zero. In the meson- 
exchange model, one finds a+ = —1.7 ■ 10~^/M7r. This small value is a consequence of 
a cancelation between the cr-exchange and the iterated p-exchange.p^ However, a priori 



there is no reason for this cancellation to remain valid when leaving the on-shell point. We 
actually observe a strong momentum dependence in the half-off-shell isoscalar T-matrix 
at threshold. The sign of this function is fully determined by the fit to the on-shell data. 

'^^Note that on the potential level, the p only contributes to the isovector channel. 
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The drawback of this approach is the blatant violation of crossing symmetry which is due 
to the iteration in the Lippmann-Schwinger equation to generate the full T-matrix. To 
be precise, the iteration generates contributions with powers in v which should be zero if 
crossing were to be respected. Some of those coefficients turn out to have non-negligible 
strength. However, one can calculate the subthreshold parameters, i.e. the expansion 
of the invariant amplitudes around the point i/ = t = in the v — t plane, which are 
allowed by crossing. For doing that, we used for convenience model 1 of ref. These 
are tabulated in table 1 (for definitions, see ref. 



Amp. 




Exp. 




Exp. 


3^02 


Exp. 




-1.54 


-1.46 ±0.10 


1.20 


1.14 ±0.02 


0.018 


0.036 ±0.003 


A-jv 


-11.46 


-8.83 ±0.10 


-0.28 


-0.374 ±0.002 


-0.028 


-0.015 ±0.002 


B+/u 


-2.6 


-3.54 ±0.06 


0.13 


0.18 ±0.01 


-0.004 


-0.01 


B- 


12.64 


10.36 ±0.10 


0.20 


0.24 ±0.01 


-0.02 


0.025 ±0.002 



Table 1: Subthreshold parameters calculated in the meson-exchange model after ne- 
glecting the crossing-violating terms (which affects in particular A~ /v and /v). The 



data are from ref. ||20| and the units are appropriate powers of the inverse pion mass. 



4. We now turn to the reaction pp ppir^. The pertinent transition matrix element is 
calculated in the distorted wave Born approximation, allowing us to properly include the 
final state nucleon-nucleon interaction. This was shown to be crucial to get the correct 



energy-dependence of the reaction under consideration |21|. 

As it was stated above, it is only the isoscalar T-matrix that enters in this reaction. 
It turns out, however, that the results of the two approaches described before are very 
different. The chiral tree calculation leads to a destructive interference with respect to 
the direct contribution (cf. fig. |I]a) (see e.g. 0), whereas in the meson-exchange model 
one finds a constructive interference (the main difference of the curves presented here 
to the ones in ||^ is that now the backward-in-time rescattering is considered in addition; 
details will be described elsewhere [0). The results for the different approaches are shown 
in fig. 1^ (using the LECs given in Eq.(^) for the chiral calculation). In the numerical 
treatment of the chiral tree level approach we deviate from the strict CHPT treatment of 
the production operator to the extent that we apply a form factor on the pion-nucleon 
vertex. This was done not only to achieve better convergence of the integral but also to 
prevent the production operator to show the abovementioned unphysical linear growth 
in the momentum. Although this treatment changes the quantitative results it does 
not effect their qualitative behavior. We have checked that while the actual strength 
of the rescattering contribution in CHPT is very sensitive to the values of the LECs 
Cj, its sign is uniquely fixed to this order (dimension two). In both results there is a 
remaining discrepancy to the data. A possible mechanism to fill this gap are the so called 
heavy meson exchanges 0. However, their sign is such that they give a constructive 
interference with the direct term (actually allowing for a good description of the data 
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|22|). Heavy meson exchange will therefore even worsen the discrepancy for the tree level 
chiral calculation 

We also applied for the first time the CHPT approach to the reaction pp — > dTT~^. Here 
neither the direct contribution nor the heavy meson exchanges contribute considerably 



23| . In this reaction a charge exchange is allowed and thus the much larger isovector 
channel of the vrA^ system can contribute. To make the comparison of the different ap- 
proaches in this reaction, we show in fig. ^ the effect of including the isoscalar rescattering 
with respect to the isovector contribution. As before, the interference pattern of the two 
approaches is totally different. While in case of the meson-exchange picture the isoscalar 
rescattering is bringing the theory in agreement with the data, using tree level chiral 
perturbation theory worsens the description. We conclude from this comparison that 
more detailed studies of this reaction might eventually lead to a deeper understanding 
of the relative size of the isovector and isoscalar components entering the A^A^ —>■ NNtt 
transition operator. 

5. To summarize, we have presented a detailed study about the relative sign of the 
rescattering and direct terms for single pion production in pp collisions in the framework 
of meson-exchange and tree level dimension two chiral perturbation theory calculations. 
While the latter is a more fundamental approach, the underlying (small) isoscalar scat- 
tering amplitude and the related A^A^ NNtt transition operator are not sufficiently 
accurate to the order they have been treated so far for the momenta involved. We can 
make this statement even more precise. Only at one loop order, terms proportional to the 
squared momentum of the exchanged pion appear in the nN amplitude. We expect exactly 
these terms to control the ttN T-matrix at large space-like momenta. Neglecting these 
terms is most probably at the heart of the sign discrepancy discussed here. We stress again 
that this does not mean that chiral perturbation theory can not be used, but rather that 
one has to work out higher order (loop) effects before one can draw decisive conclusions. 
In any case, a sufficiently improved CHPT calculation will have to be supplemented by 
heavy meson exchanges (or contact terms representing this short-range physics) to be 
able to describe the pp — ppn^ data. The more successful meson-exchange model, which 
describes the on-shell vrA^ scattering data and does much better for pp — ^ ppn^, has the 
deficiency of violating crossing symmetry. One can therefore not definitively conclude 
that this is the reason for the relative sign difference to the chiral approach, although 
this appears to be a fairly exotic possibility. In that respect, more attention should be 
focused on pp dn^ since it is much less sensitive to the heavy meson exchanges but 
still sensitive enough to the isoscalar scattering amplitude. 
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Figures 




Figure 1: Different contributions considered for the meson production, a), b) and 
c) are referred to as the direct, rescattering and heavy meson-exchange terms, in 
order. NN interactions are depicted by the blobs. 
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Figure 2: Comparison of the tvN phase shifts Su and 5*31 for the meson-exchange 
model (sohd hnes) and the tree level CHPT calculation based on parameter sets 1, 
Eq.(3) (dashed line) and 2, Eq.(4) (dot-dashed lines) with the data. 
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Figure 3: Results for S-wave pp — > pp7r°. The long-dashed hue is the result of the 
rescattering contribution from the meson exchange-model. Adding the direct con- 
tribution to this gives the upper solid line (both interfere constructively as indicated 
by the arrow). The short-dashed line is the result of the rescattering contribution 
from tree level CHPT (parameters from Eq.(3)). Adding the same amplitude for 
the direct contribution as before leads to the lower solid line (the arrow indicates 
the destructive interference). Note that heavy meson exchanges are not included. 77 
is the maximal pion momentum in units of the pion mass. 
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Figure 4: Influence of tlie isoscalar nN amplitude on S-wave pp — > dir'^. The long- 
dashed line shows the isovector rescattering of the meson-exchange model. Adding 
the isoscalar contribution leads to the upper solid line. The short-dashed hue is 
the isovector rescattering contribution of the tree level CHPT. Adding the CHPT- 
isoscalar amplitude to this leads to the lower solid curve, which substantially deviates 
from the data. Note that the effects of heavy meson exchanges and the direct term 
are small and not included here. In this case, rj = |g |/M^. 
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